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Thermal conduction in accretion disk coronae 

Andrzej Maciolek-Niedzwiecki^, Julian H. Krolik^ and Andrzej A. Zdziarski^ 

ABSTRACT 

We study the effects of thermal conduction in a hot, active corona above an accretion 
disk. We assume that all of the dissipative heating takes place in the corona. We find that 
^\ I the importance of conduction decreases with increases in the local dissipative compactness of 

0^ ' the corona, ^dissjoc, and increases with increasing abundance of electron-positron pairs. For 

•^diss.ioc < 1, a significant fraction of the energy released in the corona may be carried away by 
i-j^ I the conductive flux, leading to formation of a relatively hot transition layer below the base of 

[•-[- , ' the corona. Comptonization of disk radiation in such a layer may account for the presence of 

. _^ . soft X-ray excesses in the spectra emitted by disk-corona systems. 

^ ' Subject headings: accretion, accretion disks — conduction — galaxies: active — X-rays: galaxies 
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An accretion disk with a dissipative, thermal corona has recently become a popular model for 
interpreting the UV/X/7-ray spectra from Seyfert galaxies and Galactic black hole candidates (see, e.g., 
Svensson 1996). Recent studies have focused on the correct description of the spectrum from such a system 
(e.g.. Stern et al. 1995, Poutanen & Svensson 1996, Sincell & Krolik 1997) whereas relatively little effort 
has been taken to achieve a physically self-consistent internal structure of the corona. In particular, most 
^ I of the existing models of disk-corona systems involve two simplifying assumptions: isothermal structure of 

C^ ' the corona, and purely radiative energy exchange between the corona and the disk. We relax the above 

assumptions by considering effects of thermal conduction, which will obviously occur in a hot corona 
K^ I situated above the surface of a cold disk. Such a treatment allows us to determine self-consistently the 

; I ' electron density and temperature distributions in the corona. 
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We consider two limiting cases: a corona consisting purely of electron-positron pairs, and a corona 
without pairs (a pure electron-proton plasma). The calculations of thermal pair equilibria (e.g.. Stern et al. 
1995) show that both cases (e^ and ep dominated coronae) may be expected in accreting black hole systems, 
depending on both the compactness (defined below) and the temperature. We assume that ep coronae 
satisfy the condition of hydrostatic equilibrium. For e^ coronae, the hydrostatic model is inconsistent and 
large-scale coronal expansion is expected, as we discuss in section 2.5. Leaving the investigation of the 
exact structure of pair coronae for future studies, we consider here a constant pressure structure to see 
what effect conduction may have in the presence of pairs. 
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2.1. Notation and basic assumptions 

We assume that the corona has a slab geometry with a scale height h and a surface area A. We also 
assume a sandwich geometry (coronae located on both sides of the disk) so that we can consider only half 
of such a plane symmetrical system. The plasma in the corona is assumed to be optically thin (rtot <, 1, 
where rtot = o't L ne{z)dz is the total optical depth, He is the electron density, and z is the distance from 
the corona bottom) and to have a thermal distribution of electron energies (with the temperature varying 
across the corona). These assumptions are supported by observations of Seyfert 1 galaxies (e.g., Zdziarski 
et al. 1995; Gondek et al. 1996) and Galactic black hole candidates (e.g., Gicrlinski ct al. 1997), whose 
spectra are well fitted by optically thin thermal Comptonization. 

We measure the electron density and the temperature of the plasma by the following dimcnsionless 
parameters, 

e^^, (1) 

TOcC 
Tx = ne{x)(7Th, (2) 

where x = z/h. The latter parameter is the differential Thomson depth at x in the interval dx, i.e. 
Tx — dr/dx^ where dr — ne{z)aTdz. 

Following Haardt & Maraschi (1991, 1993), we assume that most of the accretion energy is dissipated 
in the corona. The dissipated power is lost due to thermal conduction and radiation losses. Half of the 
high energy radiation from the corona irradiates the cold disk and the other half escapes directly from the 
source (we neglect for simplicity the anisotropy of the scattering process in the corona and reflection of hard 
photons from the cold matter). The coronal radiation intercepted by the disk is reprocessed and reemitted 
at much lower energies, which provides soft seed photons for Comptonization in the corona. The dissipated 
power, idiss: a-iid both the corona and disk radiation luminosities, ih a-nd Lg, respectively, are characterized 
by compactness parameters: 
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[Lh and Lg correspond to the total power of radiation emitted by one side of the corona and the disk, 
while the observed luminosities of those layers are described by L\J2 and ^ Lsexp(— rtot), respectively.] 
Our main parameter is, however, the local dissipative compactness parameter (introduced by Bjornsson & 
Svensson 1992), 

Ij _ idiss,loc CTt , , 

-^disSjloC ^ '-, q" 7 \^) 

where £diss,ioc = Ldissh^/{Ah), and thus 

'''dissjoc '''diss~r- \^ } 

This parameter is crucial for determining the efficiency of radiative processes in the corona, since it involves 
the luminosity in a cubic volume with the size h, which is the characteristic local volume in the corona. 
Analogously, we define the local radiative compactnesses, 

^h.loc ^ ^~^' -^s.loc ^~ ^~a' ' ^ -^ 

We emphasize the importance of the difference between the local compactness and global compactness, 
i?giob = Lctt / (RiTieC^) , where R is the distance from the central source, in the corona geometry (see also 



Gondek et al. 1996). It is the local compactness which is the most relevant for parameterizing the character 
of the coronal thermal and pair balance, which fact is clear from the form of our equations (|38|-|40|) below. 
It is equally true for other problems, for example, for the issue of whether the electrons in the corona can 
achieve a thermal distribution function (cf. Ghisellini et al. 1993). The electron distribution function 
is determined by a competition between the electron-electron relaxation timescale (determined by the 
electron density and mean energy) and the electron cooling timescale (which depends on the radiation 
energy density). Their ratio is then a function of the optical depth and the local compactness, not the 
global compactness. Since ^loc is always smaller than i'giob, and sometimes very much smaller, this can be 
an important distinction. The cooling time becomes relatively longer as £ falls, so using £giob instead of £ioc 
can lead to a spurious difficulty in achieving an equilibrium distribution function. 

The importance of the conduction process will be measured by the fraction of the dissipated power 
transported out of the corona by the conductive flux, 

_ \qo\A 



where Qq is the heat flux through the corona bottom. 

Taking into account both the radiative energy exchange and the conductive transport of thermal 
energy, we can relate the corona and disk luminosities to the dissipated power, 

Lh ~ idiss(l ^ Eg), -^s — 0.5ih + Egidiss, (8) 

(we neglected here the initial energy of soft photons scattered in the corona as they contribute negligibly to 
the corona luminosity). Then 

7^-0.5(1 + 6,), J^^l-e,. (9) 



2.2. Heating 

We assume that all electrons in the corona are heated at a constant rate, 7. Then the heating rate per 
unit volume is, 

ff = 7ne(z), (10) 

which corresponds to a local dissipative compactness, 



where 



^dissjoc — r^tot, (11) 



r-^^. (12) 



2.3. Radiative processes 

Radiative cooling in the corona is due to unsaturated Comptonization and bremsstrahlung. The 
Compton cooling rate is 

Ccompt = max(46, 166^)neCCrTC/rad, (13) 



where C/rad is the radiation energy density, ric is the electron density, and the expressions in parentheses give 
the average relative photon energy gain per scattering in the nonrelativistic (40) and relativistic (160^) 
cases, respectively (for 0.2 < Q < 0.3, we use 809'^ — 328^ + 7.20 in order to achieve a smooth transition 
between the two regimes) . We take into account cooling of electrons on both soft photons from the disk 
and photons already scattered in the corona. As photons scattered many times in the corona may enter the 
Klein-Nishina regime, we include only photons satisfying the condition 

e<min(e-\l) (14) 

where e = hv/mcC^. Then, 

t'rad = -. , (15) 

c A 
where the fraction of the energy density that scatters in the Thomson regime is 

pi„(i/e,i) p^^^^^ 

F{e) is the energy spectrum of the hard radiation from the corona [we assume a power-law spectrum 
with the energy index given by eqs. (1) and (2) in Zdziarski et al. (1994) and an exponential cut-off at 
Ec = 1.6fcr], and Eq is the soft photons' energy. The bremsstrahlung cooling rate is 

Cbrom = nlmeC^aTlhicmiQ), (17) 

where 

. .^, r4^/2a/^-3/2(eV2 + l.78lel•«4), 6 < 1, ,^.. 

'^■^°'"^^^l4.5a^^-Mln(1.12e + 0.42) + 1.5], 6 > 1, ^ ' 

and af is the fine-structure constant (Svensson 1982). 

2.4. Thermal conduction 

The classical heat flux is expressed by 



<Zclass = -'«-l-7 (19) 

dz 
where k is the thermal conductivity. In a fully ionized hydrogen plasma, 

«: = 0exr^ (20) 

where the parameter (/)c allows for a reduction of the heat flux by magnetic fields, turbulence etc. (and is 
taken here to be 0c = 1) and 

•2\^/' fcV2 



(Spitzer 1962), where 



X=18 - 1/2 , , (21) 

\T^J mje'i In A 

lnA = 29.7 + ln[(re/10^K)n-i] (22) 



is the Coulomb logarithm. As In A depends logarithmically on electron density and temperature, we assume 
a constant value of InA = 30.7, which corresponds to x = 6 x lO^^erg s^^cm^^K^^/^. The classical 



expression for the heat flux is based on the assumption that the electron mean free path is short compared 
to the temperature height scale, r/|VT|. In the limit |VT| — > oo the heat flux is said to be "saturated", 
and Cowie & McKee (1977) suggested the form 

|<?sat| = bct^spcl, (23) 

where Cs = {kT/ ji)^''^ is the isothermal sound speed, jjl is the mean mass per particle {jjl ~ nip/ 2 for 
ep plasmas and m^ for e plasmas), p is the mass density and (/)s parametrizes the uncertainty of the 
appropriate value of gsat- Both experimental and theoretical results constrain (ps to the range (fig = 0.2-1.1 
with some preference for the value of (pg — 0.3 (Balbus & McKee 1982; Giuliani 1984). We use here 
4>s — 0.3 but test the sensitivity of the results to the value of 4>s. The saturated heat flux given by eq. ( p3| ) 
corresponds to advection of the plasma's thermal energy at a rate comparable to the sound speed. In order 
to smoothly implement the transition from classical diffusive to saturated trasport, Balbus & McKee (1982) 
defined the effective heat flux as the harmonic mean of ^sat and (/class, 



where 

Qclass 



(25) 
qs&t 

The equation of heat transfer in the vertical direction, 

~r = H ~ Ccompt — Cbrcm, (26) 

dz 
may then be written in a dimensionless form, 

where 

f{e,T,) = fdiss.iocT:, % ' max(49, 166^) - Tr, + Zbrcm(e)T2, (28) 

■^diss 

(29) 
(30) 

(31) 
and 6 =1 and 2 for electron-positron and electron-proton plasmas, respectively. 
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2.5. Hydrostatic equilibrium 

For an electron-proton corona, the condition h <ti R is satisfied and the following approximate form of 
the hydrostatic equilibrium equation may be used, 

dP GM . , , , 

^ = — ^(^ + ^d)p(-), (32) 



where P is the pressure in the corona and Zd is the distance of the corona bottom from the disk midplane. 
The pressure is the sum of radiation and gas pressures. From 



Prad 
^gas 



= 4 
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issjoc ~ 



4i.s 2Qt' 



(33) 



Qt Ki 0.1( 



U/4 



(Pietrini & Krohk 1995), and eq. (^) (with e, < 0.5, as found in section 4), we find 
that the gas pressure dominates for ^iss.ioc < 0.15. At higher compactnesses thermal conduction becomes 
unimportant (section 4), so we neglect radiation pressure. Then equation (B3) with P = 2nekT gives 



dx 



e 



e' 



1 m^ 
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(34) 



where Xd = Zd/h, r — R/R^ch a-nd i?sch = 2GM /c^ . For an isothermal corona the above equation gives, 

{h/Rf = 8r^&c, (35) 

nip 

where 6c is the (constant) temperature of the corona. In section 4 we find that thermal conduction 
significantly affects the corona structure near its base in a way that may increase the corona total optical 
depth. However, the temperature at higher altitudes is roughly constant and the density profile does not 
differ from that found in isothermal approximation. Therefore, we use the scale hight determined by eq. 
(p5|), approximating O^ by the temperature at the top of the corona, Otop (see section 2.7). 

The distance, Zd, is the sum of the disk height, /idisk, and the height of the transition layer affected by 
thermal conduction, ht\. Using the results of Svensson & Zdziarski (1994), we find that the ratio of corona 
height to disk height is h/hdisk ^ 100. The value of hti is determined by the solution of the conduction 
problem and (as we find in section 4) typically hti ^ O.lh when conduction is important. In this case, 
the value of Zd is dominated by the height of the transition layer. The exact value can be obtained with 
additional parameters of the model (accretion rate, viscosity parameter etc.) necessary for the calculation 

of /idisk- 

For electron-positron coronae, the gravitational force acting on the corona is too small to keep it 
bound because the escape speed is ~ cr~^'^, whereas the sound speed is ~ cmin(l/-\/3, 6^'^) and Compton 
equilibrium generically gives temperatures O not much less than unity. Consequently, in a pure pair plasma 
h ^ R and strong outflows can be expected unless the plasma is confined by magnetic fields. 

In order to check the importance of thermal conduction in pair coronae, we consider here a zone of 
e plasma at low altitudes above the disk (as it is unlikely that the heating mechanism operates at large 
distances from the disk) close to hydrostatic equilibrium so that the corona exists. As the pressure gradient 
induced by the pair gravity is negligible, we expect that such a zone is almost isobaric and so we will 
describe that zone by eq. (|3J) with the null gravity term (the second term in brackets). 



2.6. Transition layer 

As we find in section 4 below, the bottom of the corona is much hotter than the disk. Thus a transition 
layer between the disk and the corona must exist over which the temperature decreases from the coronal to 
disk value. Decrease of the temperature in the transition layer is accompanied by a decrease of conductive 
heat flux, as the energy transported from the disk is gradually radiated away. At some point the conductive 



flux vanishes and the temperature levels off at a value determined by the balance between radiative heating 
and cooling. We call this point the base of the transition layer. Below the base the temperature is still 
higher than the "normal" disk temperature, due to Compton heating by the coronal radiation. We do not 
investigate properties of that Compton heated layer as atomic processes (which are neglected in this paper) 
may become important there (a detailed analysis of such regions has recently been done by Rozaiiska & 
Czerny 1996 and SinceU & Krolik 1997). 

The structure of the transition layer affected by thermal conduction may be found in the framework of 
our model with H = 0. We take into account Compton heating of electrons in the transition layer. Then 
the right-hand side of the heat transfer equation (|2^) takes the form, 



/(e,r,) = 4issjocrJo.577exp(- y r,dx) Vj^ [max(4e, 166^) - (e)] 

+ isidL [niax(4e, 166^) - (eo)] | + ?brcm(e)T2, ^gg) 

where rj is given by eq. ( |l6|) and 



J Fede 



(37) 



" jFds- 

The integrals over the energy spectrum of the radiation emitted from the corona have upper limits equal 
of min(8~"'^, 1), and (eq) = 5 x 10~^ is assumed for soft radiation from the disk. (For Galactic black 
hole candidates, (eo) = 10^'^ would be a more appropriate value, but the results depend weakly on this 
parameter.) 

While for ep coronae the transition layer will obviously consist of an electron-proton plasma, the 
composition of the transition layer for e coronae is rather uncertain and depends on the mechanism of 
formation of the transition layer. If it is formed by the heating of the outer parts of the accretion disk, we 
can expect that it is composed mostly of ep with some addition of pairs inflowing from the corona. On 
the other hand, the transition layer may be formed together with the active region and then it would be 
dominated by electron-positron plasma. We will examine here the latter possibility (by assuming that both 
the corona and the transition layer consist only of e plasma) to estimate the maximum effect of thermal 
conduction, which is achieved for a pure pair plasma (see section 4). 



3. CALCULATIONAL PROCEDURE 
3.1. Differential Equations 

Eqs. ( P7| ) and (|3J) give three first-order differential equations, 

f = e', (38) 

ax 

^ = /(©' ^-)^^^^?^ + ^0''5(^ + ^d) - 0.566'' - 2.5^, (39) 

^ = -|[Q' + .9(^ + ^<i)]' (40) 

where f{Q,Tx) is given by eqs. ( P8|) and (|3q ) for the corona and the transition layer, respectively, and 

r 20c, for proton-electron plasma, , , 

\ 0, for electron-positron plasma. 



The null value of g for e coronae results from the isobaric approximation, whereas g — {h/ R)'^mp/ {4rme) 
for ep coronae. We solve the above equations for Q{x), <d'{x) and t^. 



3.2. Parameters and Boundary Conditions 

We are solving three first-order differential equations which depend on two dimensionless parameters, 
Ttot and ^diss.ioc- Thus, in addition to choosing the parameters, we must also choose three boundary 
conditions, 0top, ©top: ^^'^ '''x at the top of the corona, which we call Ttop. 

The physically appropriate choice for the temperature gradient is obvious. There should be no heat 
flux at the top of the corona, so 

= 0. (42) 



dx 



top 



The other boundary conditions are subtler to determine. In the absence of conduction, the temperature 
throughout the corona would be determined by ^diss.ioc and rtot- This is the calculation performed in 
numerous other discussions of accretion disk coronae. However, heat conduction removes some heat from 
the corona, decreasing its temperature, and we do not know how much this affects the top of the corona 
before doing the calculation. We implicitly fix ©top by setting the temperature at the bottom of the 
transition layer to be effectively zero, or at least far below the coronal temperature. The value used to 
initiate the integration is then a parameter to be iterated on until self-consistency (see the Appendix for 
details) . 

The third boundary condition, Ttop, rnay be thought of as the density at the top of the corona. Because 
we have normalized our unit of distance to the scale height, it is implicitly fixed by rtot- In the case of an 
ep plasma, the density profile of the corona is approximately Gaussian (it would be exactly Gaussian if the 
temperature were exactly independent of height), so formally the corona extends to infinity. However, we 
must choose some finite number of scale lengths at which to bound our calculation. We do so by setting 
the differential Thomson depth at the top Ttop = 0.0001. If the corona were exactly isothermal, this would 
be achieved at Xtop = In ' (Ttot/Ttop), so by starting the integration at xtop with Ttop — 0.0001, all would 
be self-consistent. To account for the deviations from isothermality induced by conduction, we iterate on 
cctop, keeping Ttop fixed at 0.0001. In the case of a pair plasma, our constant pressure assumption translates 
to a density profile which approaches closer and closer to a step-function in the limit as the temperature 
profile approaches isothermal. We therefore start the calculation at x = 1 and assume an initial value of the 
coronal pressure. Then we adjust the value of the pressure to obtain the assumed Ttot- 

The optical depth, Ttot, and the local compactness, -^dissjoc, of the corona are the main parameters 
of the model. Xd is also a free parameter, but its value can be determined for a specific disk model, as 
discussed in section 2.5. Increase of Xd results in a slight increase of e^, due to the change of the density 
profile in the corona. The solution of the transition layer structure gives an obvious constraint on the 
value of Xd, which must exceed the height of the transition layer. We assume an initial value of e^, which 
determines is/^dms and th/^dias from eq. (||). The value of F follows from eq. (pT|). 

Summarizing, our model has three free parameters, Ttot, ^diss,ioc and Xd- In calculations, we first adjust 
0top and Stop (or Pgas for e coronae) as described above. Then the value of e^ from the solution with 
adjusted 8top and Xtop is used in the next iteration. The procedure ends when the calculated value of tq 
converges. As discussed in Appendix A, once the three free parameters are specified, there is a (narrow) 



allowed range (0top, ©top + A0) (and a corresponding range of Xtop) within which acceptable solutions 
exist. The members of this solution family are distinguished by the pressure at the bottom of the transition 
zone. 

For ep coronae we computed solutions for 0.01 < ^dissjoc < 0.4, while in the pair case we found solutions 
for 0.01 < ^dissjoc < 10. The lower hmit on £diss,ioc was suggested by both the observational results [see, 
e.g., table 1 in Done & Fabian (1989) noting that the local compactness is roughly an order of magnitude 
lower than the global compactness given there] and theoretical expectations [for slab coronae ^diss.ioc < 0.01 
corresponds to L < lO'^^iEdd (see section 5) which is unlikely in AGNs]. Note that for ^diss,ioc < 0.1 the 
corona may still be pair dominated (see, e.g., figure 1 in Stern et al. 1995). The upper limits are set by the 
character of our solutions: at higher ^diss,ioc conduction becomes negligible (see below). For both e^ and ep 
coronae, the optical depth range over which we computed solutions was 0.06 < rtot < 0.3. 



4. RESULTS 

We find that heat flux saturation puts an important constraint on the conduction problem in that it 
establishes a limit on the amount of thermal energy that can be transported out of the corona. Namely, the 
conductive heat flux at the bottom cannot exceed the saturated heat flux, 

/ \ 1/2 3 

gsat = 50,5 — ^4-re3/2. (43) 

\ ^ / (TTh 

As this limit is independent of the total amount of released power, the relative importance of the thermal 
conduction process increases with decreasing compactness. The maximum fraction of the released power 
that can leave the corona in the form of conductive flux is, 



-'^diss -^diss.loc 

As discussed below, the temperature cannot differ strongly from the Compton equilibrium value, which can 
be estimated from Gr ~ Q.l{lhllsY''^ (Pietrini and Krolik 1995). Then th/f-s < 2 [eq. (|)] implies, 

■Y^— < 0.2-— y^-^ , (45) 

-'^diss ^diss,loc 

where we assumed 0, = 0.3. As 9 ^ 1 in both Seyfert galaxies and black hole binaries (e.g., Zdziarski et al. 
1997), we find that conduction has a rather minor effect (transporting less than 5 per cent of the released 
power) for ^diss.ioc > 1 and ^diss,ioc > 0.15 for electron-positron and electron-proton plasmas, respectively. 

Figures 1 and 2 show the dependence of eq on ^diss.ioc- The relation suggested by eq. (p5|), eq ex ^Jj^^ j^^., 
appears to be valid for fdiss.ioc ^ 1 and ^diss,ioc > 0.04 for e^ and ep plasmas, respectively. For lower values 
of ^diss,ioc, the dependence rolls over to approach the classical heat flux case in the limit of very small 
compactness (in numerical calculations we obtain the classical heat flux limit by setting 0s = 1000). Figure 
3 shows the dependence of e^ on rtot- For lower local compactness parameters, e^ decreases with increasing 
Ttot due to a higher radiative cooling efficiency of the corona. For higher ^diss.ioc (in the saturation regime), 
eq is only weakly dependent on Ttot ! a-n increase of the cooling efficiency with rtot is compensated by an 
increase of the saturated heat flux (due to the increase of r^;). 

Solid curves on Figures 4 and 5 show the temperature, the dimcnsionless electron density and the 
heat flux profiles in the corona (x > 0) and in the transition layer {x < 0) obtained within our thermal 
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conduction model. The dashed curves show the solutions of the model in which conduction is neglected 
(the temperature here is determined by the equilibrium between heating and cooling). 

Saturation prevents the temperature in the corona from being much lower than the Compton 
equilibrium temperature. Namely, a decrease of the temperature results in a decrease of the radiative 
cooling efficiency, which increases the amount of energy that must be taken out of the corona by the 
conductive flux. This quickly exceeds the capability of the saturated transport. The amount by which the 
corona temperature may deviate from the Compton equilibrium value depends on the local compactness. 
At values of ^diss.ioc for which conduction is important, the difference may become significant (Figures 4 
and 5, left panels). In this case, conduction leads to an observationally interesting effect of the softening 
of the spectrum, as a lower plasma temperature is achieved while the optical depth remains unchanged. 
For higher compactness parameters, only a tiny difference is possible and then most of the coronal emission 
comes from a roughly isothermal region (Figures 4 and 5, right panels). 

The transition layer formed below the base of the corona has a relatively high temperature and may 
extend over a height of a few tenths of the corona height scale, as shown by solid curves at x < on Figures 
4 and 5. The existence of such a layer was overlooked in previous studies of the disk-corona interface (e.g., 
Shimura et al. 1995) as they neglected conduction, which appears to provide the major contribution of 
energy to that region. Indeed, by comparing the energy input to the transition layer due to the maximum 
(saturated) value of conductive flux with the energy input due to heat flux of coronal radiation, we flnd that 

^Compt/ltl 4,loc(e)Ttl77 \ M / 

where Tti is the optical depth of the transition layer and (e) is given by eq. ([3 7]). From this we find that, 
e.g., for the parameters of the model shown on the right panel of Figure 4 conduction provides 10 times 
more energy than Compton heating. The weak contribution of Compton heating is due to the fact that the 
average energy of photons emitted by unsaturated Comptonization is much lower than the average electron 
energy in that plasma. Then, if a purely radiative energy exchange between two phases of a plasma occurs, 
the irradiated phase achieves a much lower temperature than the hot one. For a disk-corona system this 
would lead to a sharp decline in the temperature profile, as shown by dashed curves on Figures 4 and 5. 
The optical depth and Compton parameter, j/ti = / (40 -I- 16Q^)Txdx, of the transition layer decrease 
with increasing local compactness (because a lower yti is needed at higher luminosity for the loss of the 
transported energy), e.g. rti > 0.1 and yt\ > 0.3 for parameters yielding Cq > 0.2, while Tti < 0.01 and 
yti < 0.05 for e, < 0.05. 

A few solutions corresponding to different 0top are shown on the left panel of Figure 5. The solutions 
belonging to the same family may have yti differing by a few percent depending on the relative importance 
of Comptonization and bremsstrahlung in the loss of transported energy (see Appendix A). Each of 
these solutions has slightly different pressure at the base of the transition layer, which corresponds to a 
dimensionless ionization parameter 

(47) 



Fh 



rP 

^-^ gas 

The actual solution of the corona structure is determined by the disk solution, in particular by the resulting 
value of the ionization parameter at the top of the disk. The value of the ionization parameter at the base of 
transition layer is relatively sensitive to the value of Otop, as e.g. a relative change of the top temperature, 
AQ/Q ~ 10^"*, results in a change of ionization parameter, AS/S ^ 10^^, while the transported power 
varies then by only Ae^/cg ~ 10^^. 
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As discussed in section 2.4, the value of (t>s is not exactly determined. Therefore, we checked the 
sensitivity of the results to the value of 0s (see Figure 6). It appears that the level of that sensitivity is 
determined by the optical depth of the corona. For higher optical depths (r ~ 0.3), the solution depends 
very weakly on the exact value of (j)g. However, for low optical depth (r < 0.1), e^ increases significantly 
with increasing (ps for (ps < 0.6, but levels off at higher 0^. 



5. DISCUSSION AND SUMMARY 

We have developed a model of thermal conduction in a hot, active region situated above the surface of 
much colder matter. The model is applicable to regions with approximately plane parallel geometry so that 
the heat conduction equation may be reduced to one dimension. It applies also to localized active regions 
with, e.g., pillbox geometry (Stern et al. 1995), except for eq. (H) in which the deficit of soft photons must 
be taken into account. 

The importance of thermal conduction in accretion disk coronae depends strongly on the local 
compactness parameter. For ep plasma the conduction has a negligible effect for ^diss.ioc > 0.1 independent 
of other parameters. Pairs increase the importance of conduction by raising the sound speed, and therefore 
the maximum (saturated) heat flux. Our approach to pair coronae is highly simplified (a pure e corona 
and an e transition layer, no energy losses in a coronal wind) and is aimed to estimate the maximum effect 
of thermal conduction. Even in this highly idealized situation, conduction is important only for ^diss.ioc < 1- 

The value of £diss,ioc, crucial for the conduction problem, depends on the luminosity and geometrical 
parameters of the active region. E.g., for an electron-proton corona extending over the total surface of the 
accretion disk {A ~ ttR^) we find from the definition of the local compactness, using eq. (|35| ) and r ~ 20, 

^diss,ioc = 54ey2J:_. (48) 

^Edd 

Then for 10~^ < L/LEdd < 1 and 0.2 < O < 1 we can expect local compactnesses ranging from 0.02 to 54. 
In this case conduction is important for L < 10~^LEdd, a condition which may apply in at least some AGN. 
On the other hand, a luminosity of ^ lO^^iEdd is typical, e.g., for the hard state of the Galactic black hole 
candidate Cyg X-1 (e.g., Phlips et al. 1996). Thus, conduction is likely to be important in at least some 
AGNs and Galactic black hole candidates. 

When conduction becomes important, a nonuniform temperature distribution is established in the 
corona. In addition, a relatively hot transition layer between the corona and the disk is formed. In this 
case, the coronal radiation has a spectrum steeper than in models without conduction (with the same rtot 
and ^diss,ioc)- Moreover, the spectrum will have a more complicated form than a simple power-law. In 
particular, we expect a soft power-law-like emission from the transition layer due to Comptonization of soft 
photons from the disk in addition to the coronal emission. The spectral index of this emission should be 
larger than the spectral index of the coronal emission, as the Compton parameter of the transition layer is 
lower than the Compton parameter of the corona, except for very low local compactness ^diss,ioc <, 0.01. The 
superposition of the transition layer component on the spectrum of the coronal radiation (the latter extends 
to higher energies due to higher temperature in the corona) is a possible explanation of the formation of 
soft X-ray excesses, observed in many AGNs at energies below ~ 1 keV (e.g., Wilkes & Elvis 1987; Turner 
& Pounds 1989). When conduction is not important, the corona is isothermal and a rather sharp transition 
occurs at the corona base to the region in radiative equilibrium. 
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A. Otop 

We discuss here the criteria for adjusting Otop- 

©top determines the fraction of energy released in the corona that must be transported out of corona 
by conduction. For a lower 0topj Comptonization is less efficient and a higher flux must be carried down at 
lower plasma conductivity (ex 0^'^), which requires higher temperature gradient. Thus the decrease of 0top 
implies both the increase of the heat flux at the corona bottom, go: and the decrease of the temperature at 
the corona bottom, &q. The total energy transported by conductive flux is radiated away in the transition 
layer. At the base of the transition layer the solution should have a vanishing conductive flux (as no flux can 
be transported further down due to very small plasma conductivity) and small temperature (determined 
by radiative equilibrium). Below we present analytic solutions for the transition layer structure in certain 
approximations to show qualitative properties of the solution, which constrain the possible value of 0top- 

We consider the classical heat flux limit approximation with constant pressure, Pgas =const, and we 
neglect Compton heating. Below we use the dimensionless heat flux and pressure, 



h^ = ae^/^e', (Ai) 



Q — Qclass' '^ Q 

P = P^^sh^ = 2eT,. (A2) 

For most of the transition layer, the temperature is sufficiently high (and the electron density sufficiently 
low) for cooling to be completly dominated by Comptonization. We consider then the heat transfer equation 
with the cooling term determined by Compton cooling, 

^ = 4er,4iss,ioc(4 + ih)/idiss- (A3) 

ax 



Eqs. (Al) and (A3), with the boundary conditions determined by the solution at the corona bottom, qo and 




00, have the following solution 

^=(^-^^'+'^^ + ^ol ' 9 = 2£Px + go, (A4) 

where £ = £diss,ioc(^s + £h)/£diss- We see that for a given value of go there is a unique value of 
Qmax _ [7gQ/(8a£P)]^/^ for which vanishing of q at the bottom of the transition layer is accompanied by 8 
decreasing to a very small value. For higher Go temperature does not decrease and if we tried to prolong 
the solution further down the temperature would increase implying a negative temperature gradient. The 
physical meaning of that class of solutions [with 8o > ©jf '^(go)] is that conductive transport from the 
disk is required, in addition to go, in order to provide sufficient amount of energy to the transition layer to 
balance the Compton cooling rate determined by Oo- On the other hand, for Oo < 8™'''' the temperature 
vanishes before the transported thermal energy can be radiated away due to Comptonization (in this case, 
the decrease of 8 to zero is accompanied with divergence of 8' fast enough to keep the heat flux at a flnite 
value). The decrease of 8 implies, however, the increase of t^ and thus bremsstrahlung becomes efficient 
as a cooling mechanism. To check whether bremsstrahlung may provide sufficiently high cooling rate for 
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the loss of the rest of the transported thermal energy we consider the heat conduction equation with the 
cooling term determined by bremsstrahlung. We use the bremsstrahlung cooling rate for the nonrelativistic 
case and neglect the second term in parentheses in eq. (|l8|), namely 1.7819^-*'', for calculational simphcity. 
Then, 

g=4rJeV^ (A5) 



where £b = 4V2a/7r ^/^. Eqs. (Al) and ( [A5| ) have the solution 



'"^(e^-e^' f'^" 



Aa ' "' \a 



1/2 

(A6) 



where qb and 8b are the h eat flu x and the temperature at the point where bremsstrahlung starts to 
dominate. For Qt, < Iq^j \/ IbP'^a,^ the energy transported by the conductive flux cannot be radiated away. 
This means that certain 0™'" exists such that for lower 8o no acceptable solution exists, because qb and 
06 are positively correlated with q^ and Sq, respectively. 

The solution of eqs. (|3q-[l0|) has the same qualitative properties as the approximation we investigated 
above. A certain range of the values of Oo exists then, (8™'", 9™'''^), for which we obtain physically 
acceptable solutions. For lower 9o, the energy that enters the transition layer cannot be radiated away and 
as a result the temperature in the layer would increase leading to the increase of Og. On the other hand, 
for higher 9o the energy delivered to the layer is too low to maintain its temperature and the temperature 
would decrease. Solutions corresponding to different 8™'" < 8o < 8™'*'' differ with respect to the relative 
importance of bremsstrahlung and Compton cooling in the loss of transported energy. 

As discussed above, the decrease of 8top implies a decrease of 8o and an increase of go- As the value 
of go is constrained by the saturation effect (see section 4), the above constraints on 8o imply the existence 
of a certain range of the values of 8top for which the correct behaviour of the solution in the transition 
layer can be obtained. In our calculation procedure (section 2.7), we determine this range by looking at the 
behaviour of the solution at the bottom end for different values of 8top- 
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Fig. 1.--- The dependence of e^ on i?diss,ioc in electron-positron plasmas for the value of ^^ = 0.3 (solid 
curves) compared to the classical heat flux approximation (dotted curve). For solid curves Ttot = 0.06, 0.1 
and 0.2 from top to bottom and for the dotted curve Ttot = 0.1. For €diss,ioc > 1, e^ ot £~i^^^ j^^ for the solid 
curves. 

Fig. 2. — The same as in Figure 1 but for electron-proton plasmas with Ttot = 0.1; e^ ex ^^j^g iq^, for 
4iss,loc > 0.04. 

Fig. 3. — eq as a function of the optical depth of the corona. 

Fig. 4. — The corona and transition layer structure obtained in our model (solid curves) and the solution of 
the model neglecting conduction (dotted curves) for electron-proton coronae with i'dissjoc = 0.01 (left panel) 
and 0.03 (right panel) and Ttot =0.1. The curves show (a) the temperature, (b) the electron density in the 
units of Tx — ne{x)aTh, and (c) the heat flux absolute value \q\ with respect to the power dissipated per 
surface area, L^iss/A. For ^diss.ioc = 0.01 solution is obtained for Otop = 0.52466. The solution gives e^ — 0.5, 
the ionization parameter at the base of the transition layer, S = 0.004, and the optical depth and Compton 
parameter in the transition layer, Tti — 0.4 and yti — 0.6, respectively. The absolute range of possible Btop 
is AG = 1.5 X 10-5, which yields Ae, = 2 x 10"'"' and AS = 1.5 x 10"". For ii^c = 0.03, the solution, with 
Gtop = 0.618353, gives €q = 0.2, S = 0.04, m = 0.2, yti ^ 0.3, AG = 1 x 10-^ Ae, = 3 x lO^^, AS == 10--\ 
The break of |g|A/Ldiss at the interface between the transition layer and the corona (x = 0) is due to the 
fact that the corona is efficiently heated by dissipation, whereas only the much smaller Compton heating 
takes place in the transition layer. 

Fig. 5. — The same as in Figure 4 but for electron-positron plasmas, ^diss,ioc = 0.1 (left panel) and 1 (right 
panel) and t = 0.1. For i?diss,ioc = 0.1, the solution gives 8top = 0.63155, eg = 0.2, S — 0.3, Tti — 0.12 and 
2/ti = 0.3, AG = 5 X 10-^ Ae, = 3 X 10"^, AS = lO^^. For inssAoc = 1, we get Otop = 0.650531, e, = 0.05, 
S = 3.6, Tti = 0.02 and ya = 0.05, AQ = 10"^, Ae, = 3 x 10"'', AS = 1.5 x IQ-'^. The three solid curves 
on the left panel show solutions corresponding to three different values of Gtop- 

Fig. 6. — The dependence of e^ on (j)s. Left panel: electron-positron plasmas, ^diss,ioc = 0.12, 0.12, 1.2, 1.2 
and Ttot = 0.06, 0.3, 0.06, 0.3 from top to bottom. Right-panel: electron-proton plasmas, ^diss,ioc =0.01, 
0.01, 0.08, 0.08 and Ttot = 0.1, 0.32, 0.1, 0.32 from top to bottom. 
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